Electronic structure and string tension of stripes in high-Tc 

superconductors 

Marcus Fleck, Eva Pavarini, and Ole Krogh Andersen 
Max- Planck- Institut fur Festkorperforschung, 
Heisenbergstrasse 1, D-70569 Stuttgart, Federal Republic of Germany 

(February 1, 2008) 

Abstract 

We present a simple, material-specific model Hamiltonian for stripes and use 

it to explain the observed differences between the angle-resolved photoemis- 

sion spectra of La2-aSr^Cu04 and Bi2Sr2CaCu208+a;5 e.g, along the nodal 

line and around (37r/4, 0). Next we compute the string tension for stripes and 

find it to be smallest in the materials with the highest observed Tc max- 
PACS numbers: 74.72.-h, 74.25.Jb, 75.30.Fv, 79.60.-i 
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Information about the electronic structure of the high-temperature superconducting 
cuprates (HTSC) is needed for a detailed understanding of their thermodynamic and trans- 
port properties, and of the superconductivity itself. Elastic neutron scattering experiments 
on La2-o.4-i/8Ndo.4Sri/8Cu04 have demonstrated that the Cu02 layers are unstable against 
the formation of one- dimensional (ID) "charge stripes" which act as domain walls separat- 
ing vr-shifted antiferromagnetic regions |l| . These findings have triggered an intense research 
on models in which the spin and charge degrees of freedom of the electrons are viewed as 
inhomogeneous in real space 

The electronic structure of ID domain walls in the two-dimensional (2D) Cu02 layers of 
HTSCs has recently attracted a lot of interest in connection with angle resolved photoemis- 
sion spectroscopy (ARPES) on La2_^Sr^Cu04 (LSCO) §, La2_^_j,Ndj,Sr^Cu04 (Nd-LSCO) 
PJ^, and Bi2Sr2CaCu2 08+a; (BSCCO) ||8|-[TT|. A characteristic feature of the measured Fermi 
surfaces (FSs) in the 2D Brillouin zone (BZ) is that they possess straight segments running 
parallel to the (0,0)-(7r,0) and (0,0)-(0,7r) directions P,pJ|. Such fingerprints of a ID band 
structure can be explained successfully within the half- filled stripe model [T^-|TB[. Accord- 
ing to this model, the FS is a straight line, \ky\ = j, in the case of stripes running in the 
y-direction. The stripe direction, x- or y-, alternates from one Cu02-layer to the next and 
between domains within the same layer. However, the ARPES observation of a strong depen- 
dence on the momentum perpendicular to the stripe, in particular the 0.2 eV-dispersion 
along fcy ^ for I < /c^ < TT not only in Nd-LSCO but also in underdoped LSCO 

PJ^, and in BSCCO even for < < n, poses a serious problem for the model. In 
addition, the observations of nodal quasi-particles (with momentum along the diagonal of 
the 2D BZ) in optimally doped LSCO and BSCCO but not in underdoped LSCO 
P and Nd-LSCO 0, are in direct conflict with the simple stripe picture Thus, ID 

confinement of spectral weight appears not to be a generic feature of the stripe phase. 

In this Letter we present the electronic structure of ID domain walls in a Cu02 layer, 
calculated with a material-specific model Hamiltonian. The resulting band structure for an 
electron moving in an effective stripe-potential provides a simple explanation for the observed 



differences in the ARPES data between Nd-LSCO and BSCCO pO| , [TT|] . In addition, we shall 
calculate kink-excitation energies and see that they tend to zero for materials with higher 
observed Tc at optimal doping. 

We start from the 2D Hubbard model: 



k,(7 R 



where £k is the bare conduction band, which is periodic with primitive translation vectors 
Gi = (27r, 0) and G2 = (0, 2-7?) , and Oko- [a-Ra] annihilates an electron with momentum k 
[at site R = (/,m)] and spin a. This model, with = —2t{coskx + cos ky) and U = 12t, 
has recently been solved numerically using the dynamical mean-field (DMF) approximation 



and a supercell approach [|Tj]. For doping levels 0.05 < x < 0.2, y-oriented half-filled stripes 
were found to be stable. Such a band with only nearest-neighbor hopping (t) is a reasonable 
approximation for LSCO, but for materials with higher Tcmax such as BSCCO, further 
hopping integrals {f, t", ...) are needed to describe both the FSs observed with ARPES in 
overdoped samples and also the energy bands calculated with the local density approximation 
(LDA) p^JT8|] . In the present letter, we study the effect of using such material-specific band 
shapes. The self-energy in the DMF approximation, Sr,(j {uj) , defined in terms of the one- 
electron Green function by: 

G^RRV (^) = + /i - Sr. (u)] 6rr- - Y: eke''^-(^-^'\ 

k 



was found in Ref. [0] to have a spin, site, and energy dependence which, for low energies 
and for the doping levels x = 1/ (4A^) = 1/8 and 1/12, we now model by: 

Sro- (uj) = ~2^^ {Q ■ H) — ujX — i6, \u!\ < 1.5t, 

Q = ^{'^-^^^)=i^N-l)g, + Ng,. (1) 

This yields a low-energy electronic structure which is the bare band, £k, subject to a magnetic 
stripe potential of width V and a mass renormalization 1 + A. In the t-J model, the latter 
is due to the formation of spin polarons. From the DMF calculations we obtain: V = 3t 



and A = 2, for x = 1/8. The spin-periodicity of the stripe potential is AN so that the 
primitive translation vectors are Ti = (4A^, 0) and T2 = (2A^, 1) in real space, and gi = 
TT (^2^, —1^ and g2 = vr (0, 2) in reciprocal space. The observed |1[ magnetic (charge) Bragg 
peak is Q (— 2Q) . The values of sin [Q • (/,0)] for / running from to respectively 7 and 
11 are: 0, ^2/2, -1, ^2/2, 0, -^2/2, 1, -^2/2 and 0, 1/2, -^3/2, 1, ... . In the limit 
when V ^ t, the approximate form (|l]) of the self-energy yields a ID band with energy 
~ [± (y/2) sin {Q ■ (Z, 0)} — /i] / (1 + A) for each row of atoms. Here and in the following, 
the upper (lower) sign is for a spin-up (down) electron. The zero-energy, spin-degenerate, 
metalhc stripe bands centered at /=0 and 2N have the dispersions e(k^^ky)={iT/2,ky)/ (1 + A) 
and e(^^/2,ky+7T)/ (1 + A) , respectively. 

Since for the interpretation of ARPES spectra we need the |k,(j) -project ion and since 
the stripe-potential takes the simple form ([l|), which couples |k,(T) to just |k — Q,a) and 
|k + Q,o"), it is convenient to transform to the |k + nQ, (T)-representation, with the integer 
n taking AN consecutive values. As a result: 

(k+nQ,(T \G {uj)\ k+n'Q,a') = 

5,,. [{uj (1 + A) + /i + i5} 1 - H],]-^ , 

where the un-renormalized Hamiltonian matrix, 

£^k-(2Af-l)Q .0 0. ±-j 

£k-Q ±-^0.0 
±T ±? . 

±\ £k+Q . 

.0 0. £k+2^Q 

is tri-diagonal and periodic. It is isomorphic with the Hamiltonian for a 4A^-atomic ring 
with on-site energies ^k+nQ and nearest-neighbor hopping integrals =F^/4. 
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For the bare band we shall use the LDA Cu-0 antibonding pda-hand as obtained by 



integrating out the high-energy degrees of freedom |T7|,p!8 



'd . (2) 



1 — 2rak 

Here, = |(cos fc^^ + cos A;^), fk = |(cos fc^^ — cos A;^), 1/c/ is the band- width parameter, and r 
andp are band-shape parameters. Whereas l/d = 2tp^/ [fi — {Sp + ea) /2] is always about 1.6 
eV, r contains the essential materials dependence and was recently observed to increase with 
Tcmax [in]- It describes the conduction-band content of the axial-orbital, a hybrid between 
Cu 4s, Cu 3(i322_i, apical oxygen 2pz, and farther orbitals. The axial orbital lowers the energy 
of the conduction band near (vr, 0) , but not along the nodal line. The parameter p describes 
buckling-induced plane-oxygen 2p2-content and is non-zero in materials such as YBa2Cu307 
and BSCCO. For a flat layer, p = 0, whereas for a layer so buckled that p > (1 — 2r)^ , 
the saddlepoint is bifurcated away from (vr, 0) [|17|. If we expand 1/(1 — 2ruk) in (0), ^k 
takes the form of a one-band Hamiltonian with t = [1 — p + o (r)] /Ad, t' = [r + o (r)] /Ad, 
t" = t'/2 + o (r) , etc.. Hence, t/ (1 + A) ~ (1 — p) x 125 meV is at the order of the exchange 
coupling J, and r gives the range of the hopping. For LSCO, r ~ 0.1 and p = 0, while for 



BSCCO and YBCO, r ~ 0.3 and p ~ 0.17 |T3,|T8 



In Fig. 0we show for 1/8 hole-doped LSCO and BSCCO the calculated bands projected 
onto |k) in the large BZ. The spectrum for LSCO agrees well with that obtained from 



the full DMF calculation with only nearest-neighbor hopping |]l4[: The two metallic stripe 
bands are separated by gaps from the below-lying three magnetic valence bands and from 
the above-lying three magnetic conduction bands. In the |k)-projection, the stripe band 
centered at /=0 and with bottom along FX is seen near X, where it is ~75 meV below the 
chemical potential fi. Proceeding now in the y-direction, this stripe band disperses upwards 
and passes /i near (vr, j). The top of the band is seen again near Y, where it is ~75 meV above 
fi. From Y towards F, the band initially stays fiat and as it starts to disperse downwards 
through /i, it looses its |k)-character. Finally, from F along the nodal direction towards 
M, the stripe band picks up |k)-character only when, near (f , f ) , it is above /i; the nodal 
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direction therefore appears to be gapped. 

In BSCCO the stripe bands indirectly overlap the magnetic valence and conduction 
bands. This is due to the pushing down by the axial orbital of the band near X and Y. 
Along the nodal line, however, the band is completely -and the ^k-Q and ^k+Q bands are 
nearly- unaffected so that the electronic structure is like in LSCO, but /i lies lower because 
it is dragged down with the bands near X and Y. As a consequence, /i straddles off the top of 
the magnetic valence band along the nodal line. Hence, if the striped phase of BSCCO exists 
in the superconducting state, it will have nodal quasiparticles contributed by the magnetic 
valence band. Near X and Y, /i lies higher with respect to the bands than in LSCO. It even 
touches the upper part of the stripe band at Y, as well as the equivalent part of the other 
(/=4-centered) stripe band at (^|7r, 0^ . This behavior agrees quantitatively with the ARPES 
measurements on BSCCO ||10| , p!T| . In particular, our calculation identifies both the extra 
electronic component observed near (^|, |^ , which arises from the antiferromagnetic regions 
separated by ID domain walls, and also the piece of FS observed near (^jTc, 0^ , which is part 
of the Z=4-centered stripe band. 

The tri-diagonal form of H]^, together with the smallness of V/8t, allows us to discuss 
these features analytically: As seen in Fig. |I], the shape of the ^k band is left fairly intact, 
except where it is crossed by the ^k-q or the ^k+q band, or where it is crossed by the £:k-2q 
or the £:k+2q band and the corresponding intermediate band, ^k-q or ^k+q, is close by. At 
such crossing points, H]^ may be truncated toa2x2or3x3 matrix. 

Along the nodal line FM, the band is crossed by the ^k+q and ^k-q bands, which 
are nearly parallel and have ^k+q > £^k-q- The gap surrounding fi in LSCO and lying just 

above fi in BSCCO is therefore primarily due to the interaction with the ^k-q band. The 

r 1 1/2 

(un-renormalized) energies are then: ± e'^ + (^/4)^ , where e± = (^k ± ^k-q) /2, and 
the anti-bonding state is part of the /=0-centered stripe band while the bonding state is 
part of a magnetic valence band. 

Next we consider k = [jTt, O) . Since the XM line is a mirror, £:k-2q = £k for all k along 
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the hne (^^v^^r, ky^ , and ^k-Q ~ when ky = 0. The Hamiltonian is thus hke that of an 

aba molecule: A non- bonding state, (|k - 2Q) - |k)) /y/2, remains at £:k-2Q = ^k, while the 

r 1 1/2 

energies of the bonding and anti-bonding states are e+ ± + 2(V^/4)^ . Whereas the 
bonding state is part of the magnetic valence band, the non- and anti-bonding states are 
part of the metallic stripe bands centered at respectively /=0 and 2N. 

An important scale of energies which contribute spectral weight perpendicular to the 
stripes is set by the transverse stripe fluctuations We have calculated the energies of 
various elementary kink excitations relevant for the string tension of the stripes by displacing 
the stripe potential and using periodic boundary conditions with an 8 x 6 cell. The smallest 
energies were found for the excitations shown in the top panel of Fig. ^. As shown in the 
bottom panel, the string tension is calculated to decrease with the range-parameter r, once 
r ^ 0.25, and to vanish for r ~ 0.4. Hence, the energy cost of creating a kink is compensated 
by the energy gain due to long-range hopping which allows for coherent escape of holes 
sideways to the string. Since r is strongly material dependent |T^, so is the string tension. 
Also, as may seem obvious from the two spectra in Fig. [1], the energy gained by stripe 
formation decreases with r. It has been observed that for a large number of hole-doped 
HTSCs, r calculated with the LDA correlates with the measured at optimal doping |]18| . 
Specifically, as r increases from ~0.1, to ~0.3, and to ~0.4, Tcmax increases from ~25K, to 
~60K, and to ~135K. Needless to say, the reason for this correlation is unknown, but the 
present results indicate that suppression of stripes might increase Tg. In addition, our results 
in Fig. ^ show that once r ^ 0.2, the string tension is softened by dimpling or buckling of 
the Cu02 planes as it occurs in bilayered YBa2Cu307 and BSCCO. Finally, we have found 
that the domain-wall geometry with the shorter perpendicular distance has the lower energy. 
This implies that domain-wall dynamics does weaken the potential barrier for pair-tunneling 
between charge stripes pO |. 

In Figs. Q and ^ we show how the kink distortions influence the single-particle spectra 
calculated for LSCO and BSCCO, respectively. The gaps, say at (^^]^7r, o) between the 
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two stripe bands and between the lower stripe band and the upper magnetic valence band, 
tend to get filled. Even for LSCO, the parallel excursion of neighboring domain walls shown 
in the upper left of Fig. ^ nearly fills the charge gap and gives strong spectral intensity 
in the nodal direction. Our LSCO spectra shown at the bottom of Fig. ^ are in accord 
with the experimental observation of FS features near (^^, |^ generated by dynamic stripes 
in Lai.85Sro.i5Cu04 One should, however, note the difference to BSCCO where nodal 

spectral intensity from the magnetic valence band exists in already the ordered stripe phase. 
Kinks thus influence the spectrum in the nodal direction and near /i far less for BSCCO 
than for LSCO. Consistent with our energy analysis we understand that transverse stripe 
fluctuations are weaker in LSCO due to the energy cost of loosing the charge gap along the 
nodal direction. Figs. |^ and ^ show that dominant FS features, such as the straight, k^;- 
directed stripe-band segments, the magnetic valence-band sheet near the middle of the nodal 
line in BSCCO, and the sheet near (^^^tt, 0^ from the 2A^-centered stripe band in BSCCO, 
are fairly robust. This flnding for the (^^^tt, 0^-sheet is in accord with recent experimental 



observations p!0| , pn]] and clarifles the questions raised in the experimental papers about the 
origin of such quasi-particle states. 

In conclusion, we have shown that the low-energy spectral properties of striped phases 
can be studied with a simple, material-speciflc model for electrons moving in an effective 
stripe potential. This model has allowed us to understand experimental features, such as: 
(i) the absence of spectral weight in the nodal direction for underdoped LSCO and 1/8- 
doped Nd-LSCO [Q, (ii) the presence of spectral weight in the nodal direction for BSCCO 
USHTTj, (iii) the low-energy photoemission states around (|vr,0) in BSCCO [|10|,|lT|, (iv) the 
fllling of the gap along the nodal direction in the dynamic striped phase of optimally doped 
LSCO 0, and (v) the enhancement of superconducting fluctuations due to stripe dynamics. 
Finally, we have seen that the string tension of a stripe tends to be smallest in the HTSC 
materials with the highest observed Tc max- 

We thank C. Bernhard, O. Gunnarsson, P. Horsch, O. Jepsen, and J. Zaanen for stimu- 
lating discussions. 
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FIGURES 

FIG. 1. Calculated single-particle spectra, A(k,uj) for k along the symmetry directions of 
the large BZ, of the y-oriented striped phase for 1/8 doped LSCO (top) and BSCCO (bottom). 
F = (0, 0) , X= (vr, 0) , Y= (0, vr) , M= (vr, tt) , S= (7r/2, 7r/2) . 

FIG. 2. Calculated energy cost for creating localized kinks in an 8 x 6 cell {I = —2, , 5 and 

m = 0, 5) for 1/8 hole doping as a function of the range parameter r. Circles are for flat layers 
{p = 0) and squares are for p = 0.17. As illustrated in the upper panel, the open (closed) symbols 
refer to kinks which do (not) conserve the distance between domain walls. 

FIG. 3. Calculated spectra of the y-oriented striped phase with kinks in 1/8 doped LSCO. The 
top and the lower-right (left) parts are for the kinks shown in the upper-right (left) part of Fig. ^. 
The figures of the lower panel show the spectral weight in the large BZ, integrated over an energy 
window extending from // to 25 meV below. 



FIG. 4. Spectra as in Fig. | but for BSCCO. 
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